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ABSTRACT 

Chemical  and  physical  properties  of 
nitroaromatic  and  nitramine  explosives  were  related 
to  their  molecular  structure  with  the  objective  of 
predicting  types,  sites  and  rates  of  reactivity  as  well 
as  prediction  of  transformation  intermediate  and  final 
products,  stability,  and  toxicity  potential. 
Understanding  these  relationships  is  intrinsic  to  the 
development  of  computational  chemistry  (CC) 
protocols  for  basic  research  into  new  compounds  as 
well  as  for  specific  applications  to  military  sites. 
Our  comparisons  manifested  new  trends  in 
fundamental  relationships  of  these  energetic 
compounds,  including  emergent  compounds  and  their 
reactivities  in  the  environment. 

Three  classes  of  nitro  compounds  of  military 
interest  were  compared  computationally  and  verified 
spectroscopically:  i)  nitroaromatic:  trinitrotoluene 
(TNT)  and  derivatives;  ii)  cyclic  nitramines: 
hexahydro-l,3,5-trinitro-l,3,5-triazine  (RDX)  and 
octahydro-l,3,5,7-tetrazocine  (HMX);  and  iii)  cage 
nitramine:  2,4,6,8, 10,1 2-hexanitrohexaazoiso- 

wurtzitane  (CL-20).  Hypothesis  and  methodology 
were  extended  to  include  emerging  nitro  compounds: 
dinitroanisole  (DNAN);  octonitrocubane  (ONC);  and 
tetranitrotetra-aminocubane  (TNT AC) . 

1.  INTRODUCTION 

1.1  Background 

Widespread  presence  and  potential  toxicity 
of  nitro  compounds,  their  derivatives,  and  products  of 
their  environmental  transformation  require  CC 
prediction  and  analysis  of  chemical/physical 
properties  in  relation  to  molecular  structure  as  well  as 
development  of  new  CC  protocols  applicable  to 
emerging  compounds  of  military  interest  [Kholod, 
Y.,  et  ah,  2008;  Qasim,  M.,  Kolod,  Y.  et  al.,  2007; 
Qasim,  M.,  Moore,  B.,  2007]. 


1.2  Goal 

The  hypothesis  that  the  environmental  fate 
and  effects  of  energetic  materials  can  be  predicted 
from  their  molecular  structures  was  applied  to 
emerging  compounds.  Thus,  the  theme  throughout 
this  research  comprised  defining  direct  relationships 
between  molecular  structure  and  quantum  chemical 
calculations  as  related  to  chemical  and  physical 
properties.  These  properties  dictate  chemical  and 
physical  reactivities,  as  in  the  scheme  below: 

Structure  ^  quantum  mechanical 
predictions  ^  physical  and  chemical  properties  ^ 
reactivities. 

2.  APPROACH 

Since  CC  in  combination  with  experimental 
verification  is  useful  both  in  proving  concepts  and 
ascertaining  the  chemically  feasible,  our 
methodology  involved  a  theoretical  approach  and 
spectroscopic  support.  Values  generated  via 
Advanced  Chemistry  Development  (ADC/LABS) 
were  added  as  supplemental  verification  and  for 
providing  context  to  our  research. 

Theoretical  study  involved  MOP  AC 
quantum  mechanical  and  classical  force  field 
mechanics  to  predict  most  likely  bond  lengths  and 
angles,  heat  of  formation,  steric  energy,  dipole 
moments,  solvent  accessibility  and  electrostatic 
potential  surfaces,  partial  charges  and 
HOMO/LUMO  energies.  Since  compounds  of  similar 
structures  were  compared,  MOPAC  gave  reliable 
information  quickly  with  accuracy  corresponding  to 
the  requirements  of  these  comparisons  [Qasim,  M., 
Moore,  B.  et  al.,  2007].  Correlations  between 
compound  chemical/  physical  properties  and 
molecular  structure  were  then  examined  to  discover 
possible  trends. 
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Spectroscopic  tools  included: 

i)  UV/Vis  to  measure  rates  and  reactions  of 
nucleophilic  /  base  reactions  of  nitro  compounds; 

ii)  Stopped  Flow  (SF)  to  follow  rates  of  reaction 
intermediates  to  alkaline  hydrolysis.  UV/Vis  and  SF 
also  verified  selected  CC  data  for  chemical  feasibility 
and  were  used  to  distinguish  trends  evinced  by 
reactivities  of  these  explosives. 

3.  RESULTS 

Both  tables  and  figures  generate 
observations  of  trends  pertaining  to  their  subject 
compounds.  For  example,  as  seen  in  Figs.  1  and  2, 
SF  shows  that  HMX  is  less  reactive  than  RDX,  and 
Table  1  reveals  that  this  reactivity  trend  can  also  be 
seen  from  the  HMX  larger  HOMO/LUMO  gap  as 
well  as  from  the  SF  rates  of  reaction  to  alkaline 
hydrolysis. 

Table  1.  MOP  AC  Quantum  Mechanical 

calculations 
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Table  2.  Advanced  Chemistry  Development 
(ACD/LABS)  physical  property  results 
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The  following  trends  in  properties  and  reactivities 
can  be  seen  from  Table  2. 

LogP:  RDX,  HMX,  TNT,  TNA  and  CL-20 
all  have  low  LogP  values,  indicating  aqueous 
solubility  and,  thus,  possibility  of  ground  water 
contamination.  ONC  and  TNTAC  have  high  LogP 
values,  indicating  hydrophobicity. 

Henry’s  Law:  TNT,  TNA  and  ONC  have 
low  Henry’s  Law  constants,  indicating  that  they  tend 
to  remain  in  the  gaseous  phase,  whereas  RDX,  HMX, 
TNTAC  and  CL-20  have  higher  Henry’s  Law 
constants,  indicating  easier  dissolution. 

Density:  As  predicted,  RDX,  HMX,  TNT 
and  TNA  have  lower  densities  than  the  cage 
molecules  ONC,  TNTAC  and  CL-20. 
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Three  different  comparisons  were  made 
within  each  class  of  compounds:  a)  nitramine  (RDX 
and  HMX);  b)  nitroaromatic  (TNT  and  TNA);  c) 
cage  nitro  compounds  (TNTAC,  CL-20  and  ONC). 

a)  RDX/HMX:  Based  on  CC  alone,  HMX  is 
less  reactive  than  RDX,  as  revealed  by  the  larger 
HMX  HOMO/LUMO  gap.  The  HMX  larger  size  is 
represented  by  its  higher  heat  of  formation,  larger 
molar  volume,  larger  parachor  and  greater  surface 
tension.  HMX  contains  an  additional  nitrogen  in  its 
ring  system  as  well  as  an  additional  nitro  group 
attached  to  that  nitrogen.  This  results  in  HMX 
having  higher  steric  energy  due  to  HMX  collapsing 
back  on  itself.  HMX  is  more  polar  than  RDX  as  seen 
by  its  larger  dipole  charge  and  higher  polarizability 
due  to  the  increased  number  of  nitrogens  in  the 
molecule.  Both  molecules  are  highly  similar  in  their 
LogP,  pKa,  index  of  refraction  and  density. 
Comparisons  of  RDX  and  HMX  charges  are 
represented  in  Fig.  3.  Noticeable  changes  are  seen 
between  the  negative  oxygen  atoms  of  RDX  and 
HMX. 


Figure  1.  Alkaline  hydrolysis  of  RDX  via  sodium 
hydroxide  (NaOH) 


Stopped  Flow:  40  ppm  HMX  +  various  NaOH 
at  229  nm 


Time  (seconds)  -  Log  Based 


Figure  2.  Alkaline  hydrolysis  of  HMX 
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Figure  3.  Comparison  of  charge  densities  for  RDX 
and  HMX 


b)  TNT/TNA:  TNA  is  more  sterically 
hindered  than  TNT  due  to  the  oxygen  of  the  methoxy 
group.  This  oxygen  also  gives  TNA  a  larger  dipole 
charge  and  causes  a  decrease  in  the  heat  of  formation. 
TNA  is  slightly  less  reactive  than  TNT  due  to 
interaction  between  the  hydrogen  atoms  of  the 
methoxy  group  with  the  oxygens  of  the  nitro  groups. 
It  is  also  less  reactive  with  environmental  media  due 
to  more  electrons  in  the  tt  system.  The  TNA  methoxy 
group  inductively  stabilizes  the  molecule.  The  TNA 
methoxy  group  also  renders  TNA  a  larger  molecule 
than  does  the  methyl  group  of  TNT.  This  is  evident 
in  larger  molar  volume,  parachor  and  average  mass 
data.  Spectroscopic  studies  were  conducted  using 
2,4-dinitroanisole  because  of  its  classification  as  an 
emerging  compound  of  military  interest. 
Nucleophilic  reactions  carried  out  via  alkali 
hydrolysis  at  various  concentrations  of  NaOH  were 
performed  using  TNT  and  DNAN  and  are 
represented  by  the  corresponding  UV/Vis  and  SF 
spectra.  Figs.  4  and  5.  SF  spectra  of  TNT  and  DNAN 
show  mirror-image  correspondence  of 
appearance/disappearance  of  final  products.  Further, 
SF  spectra  show  alkali  hydrolysis  rates  of  TNT 
reactions  are  faster  compared  to  DNAN.  Slower  rates 
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are  expected  in  DNAN  because  of  the  methoxy’s 
effect  to  the  ring  compared  to  the  methyl  of  TNT. 
Also,  the  additional  nitro  group  found  in  TNT 
increases  the  rate  of  the  alkali  hydrolysis  reaction. 


30ppm  2,4-DNAN  in  Water  with  1  .ON  NaOH  Reacted  Over  Time 


— 30ppm  2,4-DNAN  in  Water 

— Initial 

—After  15  Min 

After  30  Min 

— After  1  Hour 

— After  1.5  Hours 

— After  4  Hours 

— After  6  Hours 

— After  24  Hours 

200  250  300  350  400  450  500  550 

Wavelength  (nm) 


Stopped  Flow:  30  mg/L  2,4-Dinitroanisole 
in  water  +  various  [NaOH]  at  301  nm 


Stopped  Flow:  30  mg/L  2,4-Dinitroanisole 
in  water  +  various  [NaOH]  at  358  nm 


0.001  0.1  10  1000 
Time  (sec)-  Log  Based 


Figure  4.  UVA^is  and  Stopped  Flow  spectra  of  2,4- 
dinitroanisole  (DNAN)  comprise  Figure  4. 


30ppm  TNT  in  Water  with  1.0N  NaOH  Reacted  Over  Time 


stopped  Flow:  30mg/L  TNT  in  water  + 
various  [NaOH]  at  235  nm 


Stopped  Flow:  30  mg/L  TNT  in  water  + 
various  [NaOH]  at  450  nm 


Figure  5.  UVA^is  and  Stopped  Flow  spectra  of 
2,4,6-Trinitrotoluene  (TNT)  comprise  Figure  5 
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Figure  6.  Comparison  of  charge  densities  for  TNT 
and  TNA 

TNAN,  not  physically  available,  was  studied 
via  CC  to  compare  the  methoxy  group  on  its  ring  vs. 
the  methyl  group  on  the  ring  of  TNT.  DNAN  was 
studied  due  to  military  interest  as  to  its  possible 
replacement  of  TNT 


Charges 


C(l)  0.2563 
NCZ)  -0.1766 
CCS)  0.1673 
NC4)  -0.2147 
NC5)  0.2941 
OC6)  -0.1969 
OC7)  -0.2060 
NC8)  0.3025 
OC9)  -0.1396 
NCIO)  -0.1854 


ceil)  0.1676 
NC12)  -0.1749 
CC13)  0.1053 
OC14)  -0.1143 
NC15)  0.4012 
OC16)  -0.2228 
OC17)  -0.2115 
NC18)  0.3015 
OC19)  -0.0787 
OC20)  -0.0744 


CL-20 


C  h  a  rg  e  s 


c 

(1 

0 

.0  0  7  8 

O 

(1 

9 

1 

-0  .0  8  4  1 

N 

(2 

- 

0  .3  0  7  3 

o 

(2 

0 

-0  .4  3  1  9 

C 

(3 

0 

.0  3  8  3 

N 

(2 

1 

1 

0.5  5  0  5 

N 

(4 

) 

- 

0  .3  0  5  2 

O 

(2 

2 

1 

-0.14  2  0 

N 

(5 

) 

- 

0  .3  0  2  0 

o 

(2 

3 

1 

-0  .4  0  6  9 

C 

(6 

) 

- 

0.0105 

N 

(2 

4 

1 

0.5  6  3  3 

C 

(7 

) 

0 

.0  2  6  0 

o 

(2 

5 

-0.1059 

N 

(8 

- 

0  .2  6  5  9 

o 

(2 

6 

1 

-0  .4  4  5  2 

C 

(9 

) 

0 

.0141 

N 

(2 

7 

1 

0.5  3  2  6 

N 

(1 

0 

) 

-0  .3  0  2  3 

o 

(2 

8 

1 

-0  .0  9  6  4 

C 

(1 

1 

) 

0.0  2  6  2 

o 

(2 

9 

1 

-0  .4  3  3  4 

N 

(1 

2 

) 

0.5  3  3  1 

N 

(3 

0 

1 

-0  .3  0  4  5 

O 

(1 

3 

) 

-0.0  9  64 

H 

(3 

1 

1 

0.2816 

O 

(1 

4 

) 

-0 .4  5  2  8 

H 

(3 

2 

) 

0.2  8  0  9 

N 

(1 

5 

) 

0.5  22  4 

H 

(3 

3 

1 

0  .2  7  8  9 

O 

(1 

6 

-0.1253 

H 

(3 

4 

1 

0  .2  84  5 

O 

(1 

7 

-0 .4  1  0  6 

H 

(3 

5 

1 

0  .2  8  5  0 

N 

(1 

8 

) 

0.5  2  3  0 

H 

(3 

6 

) 

0  .2  8  0  6 

c)  TNTAC/CL-2Q/ONC:  Of  the  cage 
nitroaromatics,  TNTAC  is  more  sterically  hindered 
than  either  CL-20  or  ONC.  ONC  is  the  least  reactive, 
as  seen  by  the  HOMO/LUMO  gap.  Also,  ONC  is 
almost  impossible  to  form  as  seen  by  its  high  heat  of 
formation.  ONC  has  the  smallest  dipole  moment 
when  compared  to  the  other  cage  explosives.  This  is 
due  to  its  molecular  symmetry.  This  symmetry 
renders  ONC  susceptible  to  free  radical  attack  as 
opposed  to  electrophilic  attack.  ONC  structure 
allows  for  a  much  larger  surface  area.  However,  CL- 
20  has  greater  molar  volume  and  a  larger  parachor 
than  the  others.  TNTAC  has  much  more  surface 
tension  and  is  far  more  dense  than  the  other 
molecules.  All  three  molecules  possess  high  MM2 
and  heat  of  formation,  which  is  due  to  the  cage  effect 
and  the  crowdedness  of  the  nitro  groups  branching 
from  the  main  ring  system. 


ONC 


Charges 


cci>  o.oiw 

CC2)  -0.0011 
lNr(3)  0.4709 
0(4)  -0.3883 
0(5)  -0.0983 
C(0)  0.1094 
INr(V)  0.4257 
0(8)  -0.4041 
0(9)  -0.1139 
C(10)  0.0523 
INr(ll)  0.4029 
0(12)  -0.3202 
0(13)  -0.0745 
C(14)  0.0058 
lNr(15)  0.4823 
0(10)  -0.3900 


0(17)  -0.0992 
C(18)  0.0881 
lNr(19)  0.4284 
0(20)  -0.3909 
0(21)  -0.1108 
C(22)  0.0258 
dNr(23)  0.4542 
0(24)  -0.4012 
0(25)  -0.1218 
C(20)  0.0378 
dNr(27)  0.4025 
0(28)  -0.3109 
0(29)  -0.0853 
lNr(30)  0.4348 
0(31)  -0.1200 
0(32)  -0.3890 


Figure  7.  Comparison  of  charge  densities  for 
TNTAC,  CL-20,  and  ONC  comprise  Figure  7. 
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4.  CONCLUSIONS 
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